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We report new results for a frequently discussed gauge theory with twelve fermion flavors in the fundamental representation of 
the SU(3) color gauge group. The model, controversial with respect to its conformality, is important in non-perturbative studies 
searching for a viable composite Higgs mechanism Beyond the Standard Model (BSM). To resolve the controversy, we subject the 
model to opposite hypotheses inside and outside of the conformal window. In the first hypothesis we test chiral symmetry breaking 
(ySB) with its Goldstone spectrum, F„, the ^SB condensate, and several composite hadron states as the fermion mass is varied 
in a limited range with our best effort to control finite volume effects and extrapolation to the massless chiral limit. Supporting 
results for ^SB from the running coupling based on the force between static sources and some preliminary evidence for the finite 
temperature transition are also presented. In the second test for the alternate hypothesis we probe conformal behavior driven by 
a single anomalous mass dimension under the assumption of unbroken chiral symmetry. Our results show a very low level of 
confidence in the conformal scenario. Staggered lattice fermions with stout-suppressed taste breaking are used throughout the 
simulations. 
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1. Introduction 

New physics at the Large Hadron Collider could be discov- 
ered in the form of some new strongly-interacting gauge theory 
with a composite Higgs mechanism, an idea which was out- 
side experimental reach when it was first introduced as an at- 
tractive BSM scenario fl}[9|. The original framework has been 
expanded by new explorations of the multi-dimensional theory 
space of nearly conformal gauge theories IfTOUTTl where sys- 
tematic non-perturbative lattice studies play a very important 
role. New experimental results at the Tevatron ITSl . boldly in- 
terpreted as Technicolor [19], will further stimulate lattice ef- 
forts to provide a well-controlled theoretical framework. Inter- 
esting models require the theory to be very close to, but below, 
the conformal window, with a running coupling which is al- 
most constant over a large energy range. The non-perturbative 
knowledge of the critical flavor NJ" separating the two phases 
is essential and this has generated much interest and many new 
lattice studies |20H44l . 

We report new studies of an important and frequently dis- 
cussed gauge theory with twelve fermion flavors in the funda- 
mental representation of the SU(3) color gauge group. With 
Nf = 12 being close to the critical flavor number, the model 
has attracted a great deal of attention in the lattice community 
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and off-lattice as well. To establish the chiral properties of a 
gauge theory close to the conformal window is notoriously dif- 
ficult. If the chiral symmetry is broken, the fundamental param- 
eter F of the chiral Lagrangian has to be small in lattice units a 
to control cut-off effects. Since the chiral expansion has terms 
with powers of NfMl/l6n 2 F 2 , reaching the chiral regime with 
a large number of fermion flavors is particularly difficult. The 
range of aM„ values where leading chiral logs can be identified 
unambiguously will require simulations in very large volumes 
which are not in the scope of this study. We will make a case 
in this report that qualitatively different expectations inside and 
outside the conformal window allow tests of the two mutually 
exclusive hypotheses without reaching down to the chiral logs 
at very small pion masses. 

Below the conformal window, chiral symmetry is broken at 
zero fermion mass with a gap in the composite hadron spectrum 
except for the associated massless Goldstone multiplet. The an- 
alytic form of the chiral Lagrangian as a function of the fermion 
mass can be used to detect chiral log corrections, or to differ- 
entiate from conformal exponents in the transitional region be- 
fore the chiral logs are reached at low enough Goldstone pion 
masses. Approximations to gauge theories with^SB, like their 
effective Nambu-Jona-Lasinio description in the large N limit, 
are consistent with this analysis. In sharp contrast, the spec- 
trum inside the conformal window is gapless in all channels in 
the chiral limit and the scale dependence of physical quantities, 
like the fermion mass dependence of composite operators and 
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their correlators, is governed by the single critical exponent y. 



2. Tests of the^SB hypothesis 



The two competing hypotheses can be tested in search for 
chiral properties of gauge theories. There is a fundamental 
difference between the two hypotheses as implied by their re- 
spective spectra. ^ SB creates a fundamental scale F in the the- 
ory separated from the composite hadron scale with its residual 
baryon gap in the chiral limit. The pion mass can be varied from 
the ^SB scale F to the hadron scale with a transition from the 
chiral log regime to a regime without chiral analysis. The con- 
formal phase has no intrinsic scale. With ^SB this is expected 
to lead to fermion mass dependence of the spectrum in the chi- 
ral log regime, or above it, quite different from the conformal 
behavior which is very tightly constrained near the chiral limit 
of the spectrum with a single critical exponent y in the absence 
of any intrinsic scale. In a regime where lattice cutoff effects are 
negligible, this difference should be sufficient for tests whether 
the chiral loop expansion is reached, or not, on the low scale F. 

In Sections 2 and 3 we present new results for the gauge 
model with twelve fermions in the fundamental representation. 
A new kind of gauge dynamics is expected to appear at inter- 
mediate distances with walking gauge coupling or a conformal 
fixed point. This has remained controversial with recent ef- 
forts from five lattice groups I120H3011 . We have made consid- 
erable progress to resolve the controversies including tests of 
the chiral condensate and the spectrum which favor chiral sym- 
metry breaking with unusual chiral dynamics. Applying the 
XSB hypothesis to the Goldstone pion, F n , the chiral conden- 
sate, and the stable nucleon state collectively leads to a result 
of^- 2 /dof = 1.22 representing a high level of confidence. With 
the conformal hypothesis we find ^- 2 /dof = 8.79 representing 
a very low level of confidence. Applying a global analysis to 
all states we measured, the contrasting behavior is somewhat 
less dramatic but remains significant. New results on the run- 
ning coupling from the static force and our simulation of a rapid 
finite temperature transition in Polyakov loop distributions, re- 
ported elsewhere and expected in association with ^SB and its 
restoration, provide further support for our findings. In Section 
4 we will briefly summarize our conclusions with outlook for 
future work. 

We have used the tree-level Symanzik-improved gauge ac- 
tion for all simulations in this paper. The conventional /3 = 6/g 2 
lattice gauge coupling is defined as the overall factor in front of 
the well-known terms of the Symanzik lattice action. Its value is 
f3 = 2.2 for all simulations reported here for the Nj = 12 model. 
The link variables in the staggered fermion matrix were expo- 
nentially smeared with two stout steps [53]; the precise defini- 
tion of the action is given in |54|. The RHMC and HMC algo- 
rithms were deployed in all runs. For the molecular dynamics 
we made use of multiple time scales ll55l and the Omelyan in- 
tegrator ll56l . Our error analysis of hadron masses which com- 
bines systematic and statistical effects follows the frequentist 
histogram approach of the Budapest-Marseille-Wuppertal col- 
laboration |57j. The topological charge was monitored in the 
simulations with frequent changes observed. 



The chiral Lagrangian for the Goldstone spectrum separated 
from the massive composite scale of hadrons exhibits, order by 
order, the well-known analytic form of powers in the fermion 
mass m with non-analytic chiral log corrections generated from 
pion loops close enough to the chiral limit (58). The exact func- 
tions F K {m) and M„(m) will be approximated by an analytic 
form in powers of m which is expected to hold over a limited m 
range when the Goldstone pion is in transition from the chiral 
log regime toward the composite hadron scale. Although this 
procedure has some inherent uncertainty without the chiral logs 
directly reached in simulations, its sharp contrast with the non- 
analytic fermion mass dependence of the conformal hypothesis, 
governed by the single exponent y, is sufficient to differentiate 
the two hypotheses. 

First, we will illustrate the fitting procedure with results on 
the Goldstone spectrum, F n , and the chiral condensate. This 
will be extended to the nucleon and some other composite 
hadron channels to probe parity degeneracy and residual masses 
in the chiral limit. 

2.1. Goldstone spectrum and F n fromx^> 

Figure 1 shows the Goldstone pion and F„ as a function of 
the fermion mass m in the range where we can reach the in- 
finite volume limit with confidence. The power functions of 
the fitting procedure in m contain the analytic contributions of 
the fourth order chiral Lagrangian to M n and F n . Although we 
could fit the pion spectrum with the logarithmic term included, 
its significance remains unclear. The rapid variation of F„ with 
m clearly shows that we would need a dense set of data in the 
m = 0.003 - 0.01 range to reach chiral logs at this gauge cou- 
pling. This requires lattice volumes well beyond the largest size 
48 3 x 96 which we could deploy in our simulations. 

Efforts were made for extrapolations to the infinite volume 
limit. At the lowest three m values, for finite volume corrections 
to M n and F„, and for all other states, we used the form 



M X (L S , rf) = M„ 
F„{L S , rf) = F„ 



1 M 2 _ 
1 + 2N f l6^F 2 ■ 8liA ^ 



N f M 2 



(i) 

(2) 



where ~g\ (A, rf) describes the finite volume corrections with A = 
M n ■ L s and aspect ratio rj = L t /L s from the lightest pion wrap- 
ping around the lattice and coupled to the measured state |59|. 
The form of ~g\ (A, rf) is a complicated infinite sum which con- 
tains Bessel functions and requires numerical evaluation. Since 
we are not in the chiral log regime, the prefactor of the ~g\ (A, rf) 
function was replaced by a fitted coefficient. The leading term 
of the function^ (A, if) is a special exponential Bessel function 
K\{A) which dominates in the simulation range. The fitting pro- 
cedure could be viewed as the approximate leading treatment of 
the pion which wraps around the finite volume, whether in chi- 
ral perturbation theory, or in Liischer's non-perturbative finite 
volume analysis [60 1 which does not require the chiral limit as 
long as the pion is the lightest state dominating the corrections. 
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Table 1 : Measured masses and F n with the three largest volumes in the m = 0.01 - 0.02 range and the largest volume for m > 0.02. Asterisks indicate L, = 32 when 
different from the spatial volume of the second column. M pmlc is the mass of the nucleon's parity partner. 



mass 


lattice 


M„ 


F„ 


Mis 


M sc 


M,j 




«™ pnuc 


Mmggs 


M rh „ 


M A i 


0.0100 


48 J x 96 


0.1647(23) 


0.02474(49) 


0.1650(13) 


0.16437(95) 


0.1657(10) 


0.3066(69) 


0.3051(81) 


0.247(13) 


0.1992(28) 


0.2569(83) 


0.0100 


40 J x 80 


0.1819(28) 


0.02382(39) 


0.1842(29) 


0.1835(35) 


0.1844(44) 


0.3553(93) 


0.352(16) 


0.2143(81) 


0.2166(73) 


0.237(12) 


0.0100 


32 J x 64 


0.2195(35) 


0.02234(46) 


0.2171(31) 


0.194(10) 


0.195(11) 


0.386(16) 


0.387(22) 


0.2162(53) 


0.239(19) 


0.246(21) 


0.0150 


48 3 x 96 


0.2140(14) 


0.03153(51) 


0.2167(16) 


0.2165(17) 


0.2185(18) 


0.3902(67) 


0.3881(84) 


0.296(13) 


0.2506(33) 


0.3245(87) 


0.0150 


40 3 x 80 


0.2200(23) 


0.03167(53) 


0.2210(21) 


0.2218(30) 


0.2239(34) 


0.4095(84) 


0.411(10) 


0.291(11) 


0.2574(36) 


0.327(14) 


0.0150 


32 3 x 64 


0.2322(34) 


0.03168(64) 


0.2319(11) 


0.2318(17) 


0.2341(16) 


0.4387(60) 


0.4333(84) 


0.2847(33) 


0.2699(41) 


0.324(16) 


0.0200 


40 J x 80 


0.2615(17) 


0.03934(56) 


0.2736(22)* 


0.2651(8) 


0.2766(42)* 


0.4673(62) 


0.4699(66) 


0.330(17) 


0.3049(28) 


0.361(32) 


0.0250 


32 3 x 64 


0.3098(18) 


0.04762(53) 


0.3179(17) 


0.3183(18) 


0.3231(20) 


0.563(12) 


0.563(14) 


0.4137(88) 


0.3683(19) 


0.469(14) 


0.0275 


24 J x 48 


0.3348(29) 


0.05218(85) 


0.3430(18) 


0.3425(25) 


0.3471(26) 


0.609(21) 


0.628(23) 


0.460(16) 


0.4050(69) 


0.523(34) 


0.0300 


24 J x 48 


0.3576(15) 


0.0561(11) 


0.3578(15)* 


0.3726(29) 


0.3790(40) 


0.640(12)* 


0.633(16)* 


0.470(15) 


0.4160(26)* 


0.5222(90)* 


0.0325 


24 J x 48 


0.3699(66) 


0.0588(15) 


0.3790(34) 


0.3814(62) 


0.3879(62) 


0.680(18) 


0.686(26) 


0.500(21) 


0.4481(39) 


0.548(31) 


0.0350 


24 J x 48 


0.3927(17) 


0.06422(57) 


0.4065(18) 


0.4074(19) 


0.4149(26) 


0.703(28) 


0.741(20) 


0.538(30) 


0.4725(64) 


0.669(65) 



Goldstone pion fit linear fit m=0.01 finite volume fit 




m m L 

Figure 1 : The Goldstone pion and F n from chiral symmetry breaking are shown with the fitting procedure described in the text. A representative finite volume fit is 
also shown. The infinite volume limit of M n was used in fits to F n and other composite hadron states, like the nucleon. 

i5Pion fit plotted with pion, ijPion, and scPion ijPion fit scPion fit 




Figure 2: The non-Goldstone pion spectrum is shown. The composite left plot displays the i5Pion data and fit together with fits to the Goldstone pion (magenta), 
i5Pion (solid blue), scPion (black), and ijPion (cyan). 



The M„L S > 4 lore for volume independence is clearly not ap- 
plicable in the model. We need M„L S > 8 to reach volume in- 
dependence. The infinite volume limits of M„ and F„ for each 
m were determined self-consistently from the fitting procedure 
using Eqs. (1,2) based on a set of L s values with representative 
fit results shown in Figures 1 and 4. In the higher m range fi- 
nite volume effects were hard to detect and even for the lowest 
m values sometimes volume dependence was not detectable for 
the largest lattice sizes. 

Non-Goldstone pion spectra, quite different from those found 
in QCD, are shown in Figure 2 using standard notation. They 
are not used in our global analysis. The three states we desig- 
nate as i5Pion, ijPion and scPion do not show any noticeable 



taste breaking or residual mass in the m — > chiral limit. The 
scPion is degenerate with the i5Pion and both are somewhat 
split from the true Goldstone pion. The ijPion state is further 
split as expected but the overall taste breaking is very small 
across the four pion states. This is a fairly strong indication 
that the coupling constant f3 = 2.2 where all runs are performed 
is close to the continuum limit. A very small residual mass at 
m — is not excluded for some non-Goldstone pion states de- 
pending on the details of the fitting procedure. 

The staggered meson and baryon states and correlators we 
use are defined in [61 1. For example, what we call the scPion 
and the fo meson are identified in correlator I of Table 1 in [61 1. 
Similarly, the i5Pion is from correlator VII, the ijPion is from 
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-3 cim - c 2 m 2 plotted 



{**) = c + ci • m + c 2 • m 2 (fitted) 
c Q = 0.00282 ± 0.00021 
C 1= 13.251 ±0.023 
C 2 = -4.21 ± 0.56 



X /dof= 0.81 

m fit range: 0.01 - 0.035 



0.01 



0.02 
m 



0.03 



.x 10 



-3{**) - m • Xcon - Ci C on m - C 2 con HI 2 



0. 



(1 - m va ^)<*M/) = {**) - m • xcon = 
cocon+ Ci C on • m + c 2c on • m 2 (fitted) 

X 2 /dof= 1.16 



Ocon 
Icon 



0.00229 ± 0.00036 
-0.193 ±0.039 
59.18 ±0.98 



2con 

m fit range: 0.01 - 0.035 



-i — f — i- 



0.04 



0.01 



0.02 
m 



0.03 



0.04 



Figure 3: The chiral condensate is shown on the left. After the coefficients of the quadratic fitting function were determined, the plot shows data for {i//ifi)-cim+C2tn 2 
for better visual display of the chiral limit at m = and its error (in magenta) coming from the fitted cq constant part. The right side is the quadratic fit to 
(if/ifr) - m ■ Xcon which is defined in the text and directly measured from zero momentum sum rules and independently from functions of the inverse staggered 
fermion matrix. The fitting function is co co „ + c\ con m + C2c 0n rr? . After the coefficients of the quadratic fitting function were determined, the plot shows data (blue 
points) for (i//ifi) - m ■ Xcon ~ ci cm m — C2 CO n'n 2 for better visual display of the chiral limit at m = and its error (in cyan) coming from the fitted cq co „ constant 
part. For comparison, the left side plot is redisplayed showing consistency between the two different and independent determinations of the chiral condensate in the 
chiral limit. For any given m always the largest volume chiral condensate data is used since the finite volume analysis is not complete. We will continue extended 
systematics at the lowest two or three m values which play an important role in the analysis. 



correlator VIII, the rho and Al mesons are from correlator III 
of Table 1 . We measure the Goldstone pion in two different 
ways, with one of them defined in correlator II of Table 1 in the 
reference. The nucleon state and its parity partner are defined 
in correlator I of Table 2 in J6T). 

2.2. Chiral condensate 

The chiral condensate (ifnft) summed over all flavors has the 
spectral representation 1621 



-2m 



-2nr 



r dAp(A) 

Jo 

f 

Ju 



m 2 + A 2 
3 dA p(A) 



A 4 m 2 + A 2 



+ c\ ■ m + C3 • m 



(3) 



where the UV-divergent integral is written in a twice-subtracted 
form in the second line [63|. The UV contribution, which is 
divergent when the cutoff aT x is removed, has a linear term 
« aT 2 ■ m and there is a third-order term as m 3 which is hard to de- 
tect for small m and survives even in the free theory limit. The 
IR finite contributions to the chiral Lagrangian have a constant 
term xBF 2 , a linear term ^B 2 -m, a quadratic term m B^F^-m 2 , 
and higher order terms in addition to logarithmic corrections 
generated from chiral loops ll64l . 

We kept a constant IR term, the linear term with UV and 
IR contributions, and the quadratic IR term in our fitting pro- 
cedure of (ljnff). The quadratic fit in Figure 3 gives a small 
non-vanishing condensate in the chiral limit which is roughly 
consistent with the GMOR [65] relation (tffiff) - 12F 2 B with 
the measured low value of F and O(l) value for B which cor- 
respond to the Goldstone pion fits in Figure 1 . The deficit be- 
tween the two sides of the GMOR relation is sensitive to the 



fitting procedure and the uncertain determination of B. The 
quadratic term in the fit is a relatively small contribution and 
trying to identify chiral logs is beyond the scope of our simula- 
tion range. 

For an independent determination, we also studied the sub- 
tracted chiral condensate operator defined with the help of the 
connected part^ co „ of the chiral susceptibility^-, 



(l-m v — )(ifnff) 
dm., 



X = "I — <<A*A> =Xc, 
dm 



(ipiff) - m ■ Xa 



(4) 



" Xdis 



Xa 



dm v 



pq \m v =m 



The derivatives d/dm and d/dm v are taken at fixed gauge 
coupling B. The derivative d/dm v is defined in the partially 
quenched functional integral of (iffiff) pq with respect to the va- 
lence mass m v and the limit m v = m is taken after differentia- 
tion. The removal of the derivative term significantly reduces 
the dominant linear part of the (ifnff) condensate. We find it 
reassuring that the two independent determinations give con- 
sistent non-vanishing results in the chiral limit as clearly shown 
in Figure 3. 

It should be noted that the M„ values in the fitting range of m 
in our analysis are below the fitting range of previous Nf - 12 
work on the chiral condensate work with considerably more un- 
certainty from using the higher range 11261 . In all fits we were on 
a fine-grained lattice in the pion mass range aM„ = 0.16 - 0.39 
and rho mass range M p = 0.2 - 0.47. In contrast, the previ- 
ous study 1 26 1 which reported conformal behavior was in the 
aM K = 0.35 - 0.67 range and rho mass range M p = 0.39 - 0.77. 
Although our new results should be made even more definitive 
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M fit M fit with M also displayed M m=0.01 5 finite volume fit 

nuc nnuc nuc r 1 nuc 




Figure 4: Nucleon and its parity partner are fitted to the constant plus linear form which are the leading contributions of the chiral Lagrangian. The blue points in 
the middle plot are the replotted nucleon data from the left to show the degeneracy of the two states. The plot on the right shows a representative finite volume fit. 



Rho meson linear fit 



Rho meson linear fit and At meson 



Table 2: The chiral condensate {ip^j) and (*//if/) - m ■ Xcom defined in the text 
and directly measured from zero momentum sum rules and independently from 
functions of the inverse staggered fermion matrix, are tabulated and used in the 
fits of Figure 3. 



mass 


lattice 


W) 


(f/nff) - m ■ Xcon 


0.0100 


48 3 x 96 


0.134896(47) 


0.006305(73) 


0.0150 


48 J x 96 


0.200647(31) 


0.012685(56) 


0.0200 


40 3 x 80 


0.266151(72) 


0.022069(76) 


0.0250 


32 3 x 64 


0.33147(10) 


0.03462(12) 


0.0275 


2A i x 48 


0.36372(40) 


0.04133(59) 


0.0300 


32 j x 32 


0.396526(84) 


0.04974(13) 


0.0325 


24 3 x 48 


0.42879(33) 


0.05781(45) 


0.0350 


2A i x 48 


0.46187(27) 


0.06807(40) 




m fit range: 0.01 - 0.035 




Figure 6: Rho meson and its splitting from the Al meson are shown. On the 
right side the magenta points reproduce the data of the rho meson from the left 
together with its linear fit. The fit parameters on the right show the linear fit to 
the At meson. 



with higher accuracy and better control on the systematics, the 
evidence is quite suggestive for a small non-vanishing chiral 
condensate in the chiral limit. 

2.3. Composite hadron spectrum in the chiral limit 



M lr quadratic fit M, ,. linear fit and scPion 




m fit range: 0.01 - 0.035 



0.01 0.02 0.03 0.04 0.01 0.02 0.03 0.04 

m m 

Figure 5: The Higgs (fo) state and its splitting from the scPion state are shown. 
The linear fit on the right works well for the Higgs (/b) state with little change 
when a quadratic term is included on the left. The blue scPion data points on 
the right and the dashed magenta fit show the fit to the scPion state. The Higgs 
will become a resonance in the chiral limit, the missing disconnected part also 
contributing, so that Higgs predictions will be challenging in future work. 

It is important to investigate the chiral limit of other com- 
posite hadron states. They further test the mass splittings be- 
tween physical states as the fermion mass m is varied and the 
measured hadron masses are subjected to chiral analysis in the 
m — > limit for important residual mass gaps above the vac- 
uum after infinite volume extrapolation. Hadron masses also 



provide useful information on parity splits in several channels. 
One composite state of great interest is the Higgs particle, if 
there is a chiral condensate close to the conformal window. We 
will briefly review new results on the nucleon state with its par- 
ity partner, the isospin partner of the Higgs (fo) state, and the 
p - Ai splitting. 

The fermion mass dependence of the nucleon and its parity 
partner is shown in Figure 4 with finite volume analysis at one 
selected fermion mass m = 0.015. The same finite volume fit 
is applied as described earlier for the pion state. The leading 
chiral linear term in the fermion mass m extrapolates to a non- 
vanishing chiral limit. The parity partner is practically degen- 
erate but this is not a surprise. Already with four flavors a near 
degeneracy was reported before by the Columbia group 1 66 1 . 

Figure 5 shows the fermion mass dependence of the Higgs 
particle without including the disconnected part of the relevant 
correlator. Strictly speaking, therefore, the state is the fo meson 
with non-zero isospin. Disconnected contributions in the cor- 
relator might shift the Higgs mass, an important issue left for 
future clarifications. Both the linear and the quadratic fits are 
shown together with the non-Goldstone scPion which is split 
down from the Higgs (Jo) state. The two states would be de- 
generate in the chiral limit with unbroken symmetry. The Higgs 
(/o) state extrapolates to a nonvanishing mass in the chiral limit 
with an Mb(j )IF ratio between 10 and 15. 

Finally, Figure 6 shows the p-meson and its A; parity partner. 
Both states extrapolate to non-vanishing mass in the chiral limit. 
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Table 3: V(R) tabulated at fermion masses m = 0.010 and m = 0.015 for lattice volume 48 3 X 96, and at m = 0.020 for lattice volume 40 3 X 80. 



m\R 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


0.010 


0.20005(49) 


0.22686(84) 


0.24638(12) 


0.26000(28) 


0.27059(55) 


0.27957(82) 


0.2872(10) 


0.2933(21) 


0.2979(42) 


0.30771(31) 


0.31250(82) 


0.015 


0.20439(21) 


0.23332(35) 


0.25270(39) 


0.26737(85) 


0.2789(17) 


0.2892(28) 


0.30214(36) 


0.3129(11) 


0.3220(31) 


0.3289(12) 


0.33576(43) 


0.020 


0.20819(39) 


0.2372(16) 


0.25961(99) 


0.27727(55) 


0.29132(74) 


0.3040(11) 
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Figure 7: V(R) data and fit for m = 0.01 is plotted on the left and comparison with perturbation theory is shown in the middle plot. The right side plot shows the 
string tension measured in nucleon mass units at m = 0.01, 0.015, 0.02 and extrapolated to the chiral limit. The finite nucleon mass gap in the chiral limit implies a 
finite string tension at m = 0. 



The split remains significant for all fermion masses and in the 
chiral limit. 

2.4. String tension and running coupling from the static force 

There are several ways to define a renormalized gauge cou- 
pling, for example, the Schrodinger Functional scheme or from 
square Wilson loops. We take the renormalized coupling as 
defined via the quark-antiquark potential V(R), extracted from 
RxT Wilson loops where the time extent T can be large. From 
the potential, one defines the force F(R) and coupling a qq (R) as 



dR R 2 



a qq {R) = 



g 2 qq (R) 

4n ' 



(5) 



The coupling is defined at the scale R of the quark-antiquark 
separation, in the infinite-volume limit L — > oo. This is differ- 
ent from the scheme using square Wilson loops, where one has 
aw(R, L) and one can choose finite R with L — > oo, or finite 
L and fixed R/L ratio. In the former case, these schemes are 
related via 

a qq (R) = a w (R)[l + 0.31551^(7?) + 0(a w (R) 2 )]. (6) 

The /} function in the qq scheme is known to 3-loops. For SU(3) 
gauge theory with Nf = 12 fundamental flavors, the location of 
the infrared fixed point to 3-loop order is a* qq = 0.3714... This 
is about 50% of the scheme-independent 2-loop value of a*, 
indicating that higher order corrections beyond 3-loop might 
not be negligible. 

A range of lattice spacings, volumes and quark masses are 
studied in the running coupling project, we show results for the 
largest volume 48 3 x 96 at /3 = 2.2 and quark masses m = 0.01 
and 0.015 and for the 40 3 x 80 run at m = 0.02. To improve the 



measurement of V(R), we use different levels of APE-smearing 
to produce a correlation matrix of Wilson loops, the lowest en- 
ergy is extracted using the generalized eigenvalue method. We 
also improve the lattice force, which is naively discretized as 
F(R +1/2) = V(R + 1) - V(R). For the Symanzik gauge ac- 
tion, the improvement is a relatively small effect, for example 
the naive value R + 1/2 = 4.5 is shifted to 4.457866... 

In Figure 7 on the left we show the measured V{R) fitted to 
the form 



V(R) = V + -+o-R. 

K 



(7) 



for m = 0.01. The m = 0.015 and m = 0.02 runs are shown on 
the right of Figure 7. For all three masses, the resulting fits are 
good, with a clear signal of linear dependence and an effective 
string tension cr. The string tension decreases with the quark 
mass, its behavior in conjunction with the mass spectrum in the 
chiral limit is under investigation and the first result is shown 
in the figure. The finite nucleon mass gap in the chiral limit 
implies a finite string tension at m = 0. 



Table 4: m lluc / V^"- 



m 


cr 


* 2 Mof 


'«nuc/ Vo 1 


0.01 

0.015 

0.02 


0.002530(81) 
0.005147(109) 
0.007189(77) 


17.1/18 
43.6/17 
6.8/14 


6.34(43) 
5.50(23) 
5.54(19) 



The renormalized coupling a qq (R) is a derivative of the po- 
tential V(R) and hence more difficult to numerically measure 
via simulations. The most accurate comparison between lattice 
simulations and perturbation theory is directly of the potential 
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V(R) itself. This is naturally given by finite potential differences 

r R a qq (R') , 
V(R)-V(Ro) = C F J ^ r L dB!, (8) 
Jro R 

where Rq is some reference point where a qq (Ro) is accurately 
measured in simulations. From this starting point, the renor- 
malized coupling runs according to perturbation theory, at some 
loop order. The result is shown in the middle of Figure 7, with 
curves at 1-, 2- and 3-loop order for the potential difference. Al- 
though progress was made in studies of important finite volume 
effects, more work is needed to bring the systematics under full 
control. It is worth noting that V(R) and its slope tend to rise at 
large R-values with increasing spatial volumes. In the current 
state of the analysis the string tension and the fast running cou- 
pling are consistent with the^SB hypothesis and do not support 
the conformal one. 

2.5. Finite temperature transition 

We present some preliminary results from our more extended 
studies of the finite temperature transition. If the ground state 
of the model has ;fSB, a phase transition is expected at some 
finite temperature in the chiral limit of massless fermions. This 
phase transition is expected to restore the chiral symmetry. If 
arguments based on universality, as implemented in a model 
framework of flavor dependence in the effective <t> 4 -theory de- 
scription, were robust as advocated l67l . the transition would be 
found to be of first order. This is not entirely clear and warrants 
careful continued investigations. 

Polyakov loop scatter plot at finite temperatures 
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Figure 8: The scatter plot of the Polyakov loop in the time direction as the 
temperature is varied using lattice sizes in a sequence of lattice sizes 48 3 x N T 
at fixed m = 0.01 with N T varied in the N T = 6, 8, 10, 12, 16 range. 

On our largest lattice, at fixed m = 0.01 and (3 = 2.2, the 
temperature was varied through an Nj sequence while the scat- 
ter plot of the Polyakov loop was monitored along the euclidean 
time direction in each run. A clear sudden transition is observed 
in the Nt =6-10 region where the Polyakov loop distribution 
jumps from the origin to a scatter plot with non-vanishing real 
part. It would be more difficult to reconcile this jump, as shown 
in Figure 8, with conformal behavior in the zero temperature 
bulk phase. 



Although we have results at other gauge couplings, and a 
variety of fermion masses as the spatial volumes and the tem- 
peratures were varied, all consistent with a finite temperature 
transition, caution is necessary before firm conclusions can be 
reached. Confirming the existence of the ^SB phase transition 
will require the m — > limits of (if/i[f) and the Polyakov loop 
distribution. The chiral condensate is a good order parameter 
for the transition. The Polyakov loop, like in QCD, could detect 
deconfinement in the transition with well-known and somewhat 
problematic interpretation issues. 

3. Testing the hypothesis of conformal chiral symmetry 

The simulation results we presented for twelve fermions in 
the fundamental representation of the SU(3) color gauge group 
favor the chiral symmetry breaking hypothesis. The pion state 
is consistent with a vanishing mass in the chiral limit and easy 
to fit with a simple quadratic function of the fermion mass. The 
non-Goldstone pion spectrum shows very little taste breaking 
at jS = 2.2 and the small splittings are consistent with expec- 
tations for staggered fermions with stout smearing. The SO(4) 
degeneracies and splittings appear to follow the pattern of QCD 
although the fermion mass dependence is significantly different. 
The fundamental scale-setting parameter F of chiral symmetry 
breaking is finite in the chiral limit. 

A non-vanishing chiral condensate is found in the chiral limit 
which is in the ballpark of the GMOR relation as suggested by 
the small value of F. We find a consistent, non-vanishing chi- 
ral limit for the subtracted chiral condensate, with the domi- 
nant linear UV-contribution removed. The nucleon states, the 
Higgs (Jo) meson, the p meson and A\ meson extrapolate to 
non-vanishing masses in the chiral limit and considerable splits 
of some of the parity partner states persist at very low fermion 
masses close to the chiral limit. There seems to be an effective 
string tension indicating confinement-like behavior below the 
string-breaking scale and the running coupling has not shown 
signs of a fixed point slowdown. In addition, there seems to be 
a rapid finite temperature transition whose nature is unclear but 
hardly favors a conformal bulk phase. Our results are consistent 
with results reported in l29l but disagree with the chiral analy- 
sis of ll26l and do not support the infrared fixed point reported 
in El. 

But is it possible that we mislead ourselves with the ^SB 
interpretation? Can we interpret the results as conformal chi- 
ral symmetry? To decide this question, a fairly stringent test 
is possible. With the conformal hypothesis the mass depen- 
dence of all physical states is controlled by the anomalous di- 
mension y for small fermion masses l35ll . Each hadron and 
F„ should scale as M n * m 1/ym and F„ * m ^ m for small m 
where y m = 1 + y. For small enough m the value of y should 
be interpreted as y* at the infrared fixed point. The chiral con- 

— hi 

densate is expected to have the behavior (ifnft) ~ c ■ m + m l *y 
when m — > 0. We selected various subsets of states for a com- 
bined fit with universal critical exponent y. We also fitted all 
measured states combined. Applying the conformal hypothe- 
sis to the chiral condensate, to F, to the pion state, and to the 
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Figure 9: The Nf = 12 chiral and conformai simultaneous fits in four channels are displayed for comparison in two select cases. 



stable nucleoli state collectively yields a total y 1 = 229 for 26 
degrees of freedom with^ 2 /dof = 8.79. This indicates a very 
low level of confidence in the hypothesis. The^SB hypothesis 
gives ^ 2 /dof = 1.22 for the same set of states. This was the 
result quoted in Section 1. The chiral and conformai fits for 
two of the four fitted states with the quoted global results are 
shown in Figure 9. Applying the global analysis to all states 
we measured, the contrasting behavior is less pronounced but 
still significant. The results disfavoring the conformai hypoth- 
esis are significant. More work is needed for higher accuracy 
and full control of the systematics, yet it is worth noting that as 
the volumes are increased at the lower quark masses, the results 
for F„ and (ifnf/) will increase or remain the same; this does not 
bode well for the conformai picture. 

4. Conclusions and outlook 

We reported new results for a frequently discussed gauge the- 
ory with twelve fermion flavors in the fundamental representa- 
tion of the SU(3) color gauge group. Our results favor with a 
significant level of confidence the ^SB scenario but close to the 
conformai window several features of this gauge theory clearly 
differ from QCD. We find a large B/F ratio which is often in- 
terpreted as strong chiral condensate enhancement. This could 
be related to the fermion mass dependence of the Goldstone 
pion and the non-Goldstone pion spectra which is different from 
what is observed in QCD with staggered fermions. In compari- 



son with QCD, we also observe significantly smaller mass split- 
tings between parity partners in several hadron channels. The 
near degeneracy of parity partner states is expected to lead to 
an S-parameter quite different from what was projected when 
QCD was scaled up by the number of flavors [68 1. To gain 
more confidence in our ^SB analysis, it is important to push 
deeper into the chiral regime and closer to the continuum limit 
than the analysis reported earlier |26|. Only future work with 
high precision simulations will be able to explain whether this 
is the source of our qualitatively different findings. 

The infrared fixed point (IRFP) of the gauge coupling re- 
ported earlier [22| is also in disagreement with the picture we 
presented. It seems to be very difficult to differentiate between 
an IRFP and a slowly walking gauge coupling close to the con- 
formal window as also indicated with MCRG studies of the 
model ||28l . It would be interesting to establish in a quantita- 
tive analysis using the Schrodinger Functional method of Il22l 
and the MCRG method of [28 1 the difference in the level of 
confidence between the two mutually exclusive hypotheses of 
IRFP or walking gauge coupling. This was precisely the main 
thrust of our work in comparing the ^SB hypothesis with the 
conformai hypothesis. 

We plan the continued investigation of the running gauge 
coupling reported in Section 2 with the ambitious goal of dif- 
ferentiating between the IRFP and walking gauge coupling sce- 
narios. If a walking gauge coupling g„ can be established ap- 
proximating IRFP behavior over some extended scale, it will be 
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important to determine the related mass anomalous dimension 
y(gfy and the S-parameter as model examples for BSM appli- 
cations. 
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